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ABSTRACT 

Aims. We investigate whether the mean star formation activity of star forming galaxies from z=0 to z=0.7 in the GOODS-S field 
can be reproduced by simple evolution models of these systems. In this case, such models might be used as first order references 
for studies at higher z to decipher when and to what extent a secular evolution is sufficient to explain the star formation history in 
galaxies. 

Methods. We selected star-forming galaxies at z=0 and at z=0.7 in IR and in UV to have access to all the recent star formation. We 
focused on galaxies with a stellar mass ranging between 10'° and 10" M Q for which the results are not biased by the selections. We 
compared the data to chemical evolution models developed for spiral galaxies and originally built to reproduce the main characteristics 
of the Milky Way and nearby spirals without fine-tuning them for the present analysis. 

Results. We find a shallow decrease in the specific star formation rate (SSFR) when the stellar mass increases. The evolution of the 
SSFR characterizing both UV and IR selected galaxies from z=0 to z=0.7 is consistent with the models built to reproduce the present 
spiral galaxies. There is no need to strongly modify of the physical conditions in galaxies to explain the average evolution of their 
star formation from z=0 to z=0.7. We use the models to predict the evolution of the star formation rate and the metallicity on a wider 
range of redshift and we compare these predictions with the results of semi-analytical models 
Conclusions. . 
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1. Introduction 

A lot of recent studies have explored the relation between the 
stellar mass and the star formation rate (SFR) in galaxies at 
different redshifts. It has been shown that star formation criti- 
cally depends on galaxy mass both at low and hi g h redshift (e.g . 
Gavazzi et ail l2002t iBrinchmann et al.L 12004 iFeulner et all 
2005L a nd reference therein ). The general process called "down- 
sizing " dCowie etaU ll996) is now commonly accepted: it can be 
summarized by an early and rapid formation of massive galaxies, 
whereas low-mass systems evolve more smoothly, being still ac- 
tively forming stars at z=0. The galaxies are usually subdivided 
into active star-forming and quiescent systems within which star 
formation has been quenched. When only star-forming galax- 
ies are concerned, the specific star formation rate (star formation 



rate divided by stellar mass) seems to exhibit a flatter distribution 
as a function of the stellar mass t han when the whole population 
of galaxies is accounted for ( e.g. Daddi et all l2007t lElbaz et all 
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I2007tllglesias-Paramo et"ail 12007^ 

Physical processes that might be at the origin of or might in- 
fluence the quenching of the star for mation in massive galaxies 
have been extensiv e ly dis cussed (e.g. Dekel & Birnboim, 2006; 
Bun dv et all 120051 120061 and references therein). Conversely 
there are only very few attempts to propose realistic star for- 
mation histories to explain the trends that are observed: most of 
the studies only explore very crude scenario s like a constant sta r 
formation rate or an in s tanta neous burst dFeulner et all 12005). 
Recently, Noe ske et al.l d2007l) have proposed a more sophisti- 
cated model with an exponential star formation history whose 
parameters are mass-dependent: less massive galaxies have a 
longer e-folding time and begin their formation more recently. 
Their model is based on two parameters: the e-folding time r and 
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redshift of galaxy formation zu which are fitted to their data. The 
redshift formation zj is likely to be representative of the bulk of 
the star formation: a galaxy with a stellar mass of ~ 10 1() Mo 
would form at z/ — 1 against Z/ = 3 for more massive objects 
~ 10 U MO. 

In this paper, we follow a different approach. We start from 
a physical model built to reproduce the mean properties of local 
star-forming disk galaxies. We analyze how this model, with- 
out any modification or adjustment, can reproduce the evolu- 
tion of the star formation observed with z. Practically speak- 
ing, we use a grid of models predicting the chemical and spec- 
trophotometric evolution of spi ral galaxies. These models were 
calibrated in the Milky Way (Boissi er" & Prantzoi, 1 19991) and 
successfully reproduce the properties of nearby spiral galaxies 
(Boissier & Prantzoi, l2000h . The redshift formation z/ is taken 
to be equal to 6 for all the galaxies and represents the time when 
the first stars begin to form. 

Our work is motivated by the recent findings that the UV 
and IR luminosity functions at intermediate redshifts (up to z 
0.7 - 1) are almost dominated by normal spiral or irregular 



galax ies (Bell et al 



120051; IZheng et al. 



2005 



2007b; 



Wolf etaU 120051; [Melbourne et all 
Zamoiski et all 120071) . The role of 



major mergers does not seem to govern the star formation at 
these redshifts, one can expect those galaxy evolution models 
that assume a smooth evolution to be appropriate for the inter- 
pretation of the observations. 

In section 2 we describe the samples of sta-forming galaxies 
at z=0 and z=0.7 and the derivation of the main parameters 
useful to this study: the star formation rates (SFR) and stellar 
masses (M stal ). The selections are performed in rest-frame IR 
(60 or 15/rni) and UV (1530A) in order to be sure to pick up all 
star-forming galaxies at these redshifts. The characteristics of 
these selections are outlined in section 3. Specific star formation 
rates (SFRs divided by M stal ) are analyzed in section 4. The 
evolutionary models are presented in section 5 where we pro- 
pose a simple analytical formulation for the SFR history of each 
galaxy. In section 6 the models are compared to the data. We 
also compare our results to more sophisticated (semi-analytical) 
models and give predictions for higher z. Section 7 is devoted to 
the conclusions. 



ate z the deepest GALEX and SPITZER surveys can be used to 
build similar samples to z=0. Practically speaking, we work at 
z=0.7 to avoid K-corrections. At this redshift the GALEX near- 
ultraviolet band at 2310 A corresponds to a 1530 A rest frame, 
the UV wavelength at which galaxies have been selected at z=0 
(far-ultraviolet band of GALEX). The SPITZER/MIPS observa- 
tions at 24 /mi correspond to ~ 15 //m in the rest frame of galax- 
ies at z=0.7. Although 15 //m does not directly correspond to the 
IRAS 60//m band, this mid-infrared wavelength range has been 
intensively studied, so we will be able to derive total infrared 
luminosities Lir from these mono-wavelength data. 

Rest-frame, near-infrared (NIR) data will be needed to esti- 
mate stellar masses, whereas the UV/IR selections ensure we can 
me asure the current st ar formation: at z=0 they were estimated 
by iBuat et alj d2007bl) from 2MASS data. At z=0.7 the IRAC 
observations at 3.6 /mi correspond to the K-band rest frame. 
All the UV da t a are corrected for Galactic extinction using the 
Schlegel et al l (1998) maps and the Galactic extinction law of 



Carde lli et al. (1989). We have gathered the main characteristics 



of the samples described below in Table 1 . 



2. 1 . z=0 samples 



IBuat et all d2007bl) built two samples of galaxies selected at 60 
/mi (IR selected sample) and at 1530A (UV selected sample) 
with a very high detection rate at 1530A and 60 /mi, respec- 
tively. A short description of the samples is made in Table 1. 
The original samples were flux-limited: the limiting luminosities 
reported in the table correspond to the faintest bins o f the lumi- 
nosity functions built with these samples (cf. Fig. 3 of Bua t et all 
2007b). Here we use the mean (volume-averaged) trends found 
from these samples. The SFRs are estimated by adding the star 
formation rate measur ed from the IR and the observ ed UV emis- 
sions as preconised by Iglesias-Paramo et alj (120061) assuming a 
constant SFR over 10 8 years. The only modification that we per- 
form h ere is to use the initial mass function (IMF) of Kroupa 
d2001l) since it is commonly used in recent studies and it is more 
consistent with the mo dels used in this paper. From Starburst99 
dLeitherer et all 1 19991) we obtain: 



Throughout this article, we use the cosmological parameters 
H = 70 km s -1 Mpc~', Q. M = 0.3 and Cl A = 0.7. All magnitudes 
are quoted in the AB system. The IR luminosity Lir is defined 
over the wavelength range 8-1000 //m. The UV luminosity Luv 
is defined as vL v . The luminosities are expressed in solar units 
withL = 3.83 lO^ergs" 1 



2. Data sets, SFRs and stellar mass estimates 

Our goal is to study the evolution of star-forming galaxies from 
the local universe to intermediate redshift (z <~ 1) by compar- 
ing observations to the predictions of models of their evolution. 
Therefore we need to build galaxy samples that are representa- 
tive of the overall star formation and we must avoid quiescent 
galaxies. The best way to focus on star-forming galaxies is to se- 
lect them according to their SFR. Newly formed stars emit most 
of their light in the UV and a large fraction of their emission is 
re-processed in the IR via dust heating. Therefore we select the 
galaxy samples by using these two wavelength ranges: UV and 
IR. 

We use local samples already built from GALEX and IRAS 
data at 1530 A and 60 /mi (Bua t etaU l2007bl) . At intermedi- 



log(SFR 1R ) Me yH = log(L IR ) Lo - 9.97 
\og(SFRw) Mo yr -i = log(L uv ) io - 9.69 



(1) 

(2) 



We follow esias-Paramo et al. (2006) to estimate the total 
SFR: 



SFR = (1 - if, ■ SFR 1R + S FRi 



(3) 



where rj is the fraction of dust emission due to the heat- 
ing by old stars and not related to the recent star forma- 
tion. The introduction of this factor was found to be m anda- 
tory for making consistent star formation ind i cators dBelll 
2003; Hirashit a et all 120031; llglesias-Paramo etail l2004l2006l). 



Following llglesias-Paramo et al. (2006) we take rj = 0.3. 

The stellar m ass is calcul ated from H ma gnitudes as in 
IBuat et alj d2007bl) . We use the lBell etall d2003l) M/L H calibra- 
tions adopting a me a n B-V color of 0.6 mag (Bu at et all 12007b) 
and a IKroupa et al] d 1993b I MF. It cor r espon ds to M/ L H =0.58 
in sola r units. According to Bell et al] d2007l) . using a I Kroupal 
d2001l) IMF instead of a IKroupa et all d 1993b one would lead 
to similar stellar m asses (within 10%). We have also checked 
dBuat et all l2007bl) that very similar masses would be obtained 
using K magnitudes instead of the H ones. 
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Table 1. Description of the samples. At z=0 all the characteristics c ome from Buat et all d200 7b). the limiting Lir is estimated from 



the luminosity at 60 yum by applying a bolometric correction of 2.5 (Bua t et aH 1200 7b). The samples at z=0.7 are fully described in 
section 2.2 





z=0 


z=0.7 




IRAS f 60 > 0.6Jy 


SPITZER/MIPS f 24 > 0.024mJy 


IR selection 


665 galaxies 


280 galaxies 




Lir > 2.5 10 9 L G 


Lir > 2.2 10 10 L o 




GALEXFUV(1530A) < 17mag 


GALEXNUV(2310A) < 25.5mag 


UV selection 


656 galaxies 


247 galaxies 




L uv > 3 10 8 L o 


L uv > 2.5 1O 9 L 



2.2. z=0.7 samples 
2.2.1. IR-selected sample 

We start with th e GOODS observations of the C handra deep field 
south (CDFS) (Dic kinson & Giavaliscol |2002|) . The source ex- 
traction and flux measurem ents were performed at 24 /jm using 
DAOPHOT dStetsorJ,[l987l) in a same way as in iLe Floc'h et all 
(2005). The detection limit at 3 cr is found to be 0.016 mJy. 
We have checked the completeness of our data by comparing 
the sl o pe of the differ e ntial c ounts with models (La gache et all 
120031: ITakeuchi et all 1200 ll) and data obtained in the HDF- 
N dCharv et all 120041) . The data appear to be complete down 
to ~0.024 mJy. Th ese values fully agree with those found by 
Elba z et al] d2007l) for the same set of da ta. We adopt a se - 
cure limit of 0.024 mJy to build our sample (El baz et all 20071) . 



Cross correlations with COMBO- 17 (Wo lf et all 120041) were 



performed to get the redshift of the sources within a tolerance 
radius of 2" dLe Floc'h et all 120051) . 88% of the sources have 
a single counterpart in 2". We restrict the final sample to these 
objects with a single cou nterpart. The optical p hotometry is re- 
trieved from the EIS data (lArnouts et all 1200 ll) 

To build the sample at z=0.7 we selected galaxies in the red- 
shift bin 0.6-0.8. We put a limit in luminosity to be complete in 
the redshift range and to avoid volume corrections. This limit is 
calculated for a limiting flux of 0.024 mJy at z=0.8, and trans- 
lates to a total infrared luminosity log(LiR) = 10.34 (L ) (see 
below for the calibration of the 24 fim flux in Lir). 280 sources 
are obtained that way (the faintest object has a 24 fim flux equal 
to 0.031 mJy). All these galaxies are detected with IRAC at 3.6 
yum. 

We must estimate the total infrared luminosity Lir from 8 to 
1000 //m to measure the "obscured" SFR. The extrapolation 
from the 24 //m emission a lone to the tota l IR emission re- 
lies on loca l templates (e.g. iLe Floc'h et all 120051 : iBell et all 
I2005L 120071) : tight relations have been found between the rest 
fram e 12-15 nm and the total infrared emission of local galax- 
ies dCharv & Elbazl I200U ITakeuchi et all l2005al) . Recent stud- 
ies based on SPITZER data suggest that there is not a strong 
evolution in the IR sp ectral energy distrib utions for interme- 
diate redshift galaxies (Zhe ng et all l2007bl) . At z=0.7, 24 /im 
corre sponds to 14 y um and a calibration at 1 5 /um can be 
used. iTakeuchi et ail d2005al) and lCharv & Elbazl d200ll) propose 
such ca librations. Here we use the relation of Takeuchi et al. 
d2005al) based on IRAS and ISO data and local templates from 
iDale et all d200ll) to be con sistent with z=0 calculations. Using 
the Chary & Elbaz (2001) relation would lead to a slightly 
higher value of Lir by 0.08 + 0.01 in log units. 
1 15 galaxies are detected at 2310 A (NUV band from GALEX). 
At z=0.7 it corresponds approximately to the rest-frame FUV 
band of GALEX centered on 1530 A (cf lBuat et alll2007al) . In 
the same way, the 3.6 yum corresponds to a rest-frame K band. 



Therefore we can avoid K-corrections. For the galaxies detected 
at 24 fim, and in UV the SFR is estimated by combining the UV 
and IR emissions in a similar way to z=0 (Eq. 3). Since galax- 
ies at higher z are more active in star formation than at z=0, we 
can expect a lower contribution of old stars to dust heating and 
hence a value of 77 l ower than 0.3. To check this issue, we fol- 
lowed the method of Iglesias-Para mo et al.l (|2004) by comparing 
SFRs calculated with Eq. 3 to those deduced from the UV lumi- 
nosity alone, co rrected for dust attenuation using the recipe of 
lBuatetalJd2005l) . The SFRs were found to be consistent if 77=0; 
i.e., all the dust heating is attributed to young stars in these ob- 
jects. Therefore we adopted 77=0 to calculate the SFRs at z=0.7. 
For the galaxies not detected in the UV, the SFR was calculated 
with the IR luminosity only (again with 77=0). The contribution 
of the UV emission at the detection limit of UV=26.2 mag is 
negligible in the estimation of the SFR (see section 4). 
As for the z=0 samples , the stellar masses were estimated fol- 
lowing B ell et al.1 ((2003) who calibrate the M/L ratios as a func- 
tion of several colors including those of the SDSS. The observed 
3.6 fim band corresponds to the K band at z=0.7 and the ob- 
served R-I color from the EIS survey (Johnson-Cousin system) 
is similar to the u-g rest-frame color (fr om the SPSS). Th erefore 
we used the u-g color-M/LK relation of lBell etaT] d2003l) for the 
305 galaxies with a R-I color. 18 galaxies have no R-I color so 
for them we used the mean M/Lk obtained for the sample i.e. 
M/Lk=0.46 (solar units). The R-I color distribution is discussed 
in the next section. 



2.2.2. UV-selected sample 



GALEX dMorrissev et all 120051) observed the CDFS for 76 ks 
in both the FUV (1530 A) and the NUV (2310 A) as part of its 
deep imaging survey. The reductio n of the data is extensively 



described in Burga reTla et al.l d2007l) . Very briefly speaking, we 
used DAOPHOT to perform PSF-fitting and disentangle close 
neighbors. The completeness at a level of 80% was obtained at 
NUV= 26.2 mag. The data are available from Bur garella et al.l 
d2007l) . The cross-i dentification wit h the C OMBO- 17 sources 
is also described in Bur garella et a l. (2007). About 70% of the 
GALEX sources are identified in COMBO- 17 with a strong de- 
pendence on the NUV magnitude. Truncating at NUV=25.5 mag 
ensures us that more than 80% of the GALEX sources are iden- 
tified. We adopt this limiting magnitude in the following and we 
restrict the final sample to objects with a single counterpart in 
COMBO- 17 ( 90 % of the UV sources have more than one coun- 
terpart in COMBO- 17). As for the IR selection to avoid volume 
corrections, we put a limit in luminosity, which ensures that we 
detect all the galaxies brighter than this limit in the redshift range 
0.6-0.8. This limit is calculated for a magnitude NUV=25.5 at 
z=0.8. It corresponds to log(Luv) = 9.39 (L ) and 247 galax- 
ies are selected. 48% (119/247) of these sources are detected at 
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24 yum. For the undetected ones, we adopted an upper limit at 
0.016 mJy, which corresponds to the detection limit at 3 cr (sec- 
tion 2.2.1). The IR luminosities Lir were estimated from the 24 
/im flux in the same way as for the IR-selected sample (section 
2.2.1). The stellar masses were estimated in the same way as for 
the IR sample at z=0.7 and 230 galaxies have a R-I color. For 
the 17 remaining objects without a R-I color, we used the mean 
M/Lk obtained for the sample i.e. M/Lk=0.42 (solar units). The 
R-I color distribution will be discussed in the next section. 

The star formation rate was estimated as for the IR selection 
when UV and 24 /mi fluxes are available. We again took 77 = 0; 
this choice is also validated by the comparison between the SFRs 
calculated with Eq. 3 and 77 = and those estimated from the UV 
luminosity corrected for dust attenuation (as for the IR selection, 
section 2.2. 1). An upper limit on the SFR was calculated for the 
galaxies with only an upper limit at 24//m . This time we could 
not neglect the contribution of the IR emission which is domi- 
nant, even at the detection limit level (see next section). 



3. Comparison of the UV and IR selections 

Before discussing the star formation activity in both samples we 
analyz e which selection (I R or UV) is best suited to our analysis. 
At z=0 lBuat etaf] (l2007bl) showed that the intrinsically brightest 
galaxies are lost in a UV selection. Conversely, intrinsically faint 
galaxies are hardly detected in IR. Because our present study is 
devoted to a comparison between low and higher z samples we 
do not discuss the intrinsically faint objects against which we 
are strongly biased (cf next section). At z=0, intrinsically bright 
galaxi es are rare and the d ifferences between both selections are 
small dBuat et all l2007bl) : however, the galaxy population seen 
either in UV or in IR is known to brighten as z increases, we 
must also check the differences in the selections at z=0.7. 

3.1. The relative contribution of the IR and UV emissions to 
the bolometric luminosity of young stars 

Atz=0|Buat et al.l(l2007bl) compared the bolometric luminosities 
of the galaxies defined as Lboi = (1 - i])Lir + Luv- Briefly sum- 
marized, we showed that the UV luminosity alone is unable to 
reproduce the bolometric luminosity even for UV-selected galax- 
ies and that the combination with the IR is mandatory. Once 
both luminosities are added to calculate Lboi, a deficiency of very 
bright objects (in terms of Lboi) is observed in the UV selection 
as compared to the IR one. Given the low number of such galax- 
ies, these differences are slight as is shown in the next section, 
which is devoted to specific star formation rates. 

In Fig. [TJ we have plotted the IR to UV luminosity ratio ver- 
sus Lboi for both samples at z=0.7. It can be seen that the IR lu- 
minosity is higher than the UV one for all the galaxies detected 
at both wavelengths since Lir/Luv > F The galaxies selected in 
UV and not detected at 24 //m exhibit the lowest bolometric lu- 
minosities with Lboi ^ 2.5 1O 1O L (log(L B0 i) < 10.4). The galax- 
ies selected in IR and not detected in UV span a very wide range 
of bolometric luminosity and hence of SFRs. The average trends 
found at z=0 are also reported in Fig.Q] The samples at z=0 and 
z=0.7 cover the same range of luminosity. The slight difference 
found in the distributions of the IR to UV l uminosity ratio has 
been disc ussed in terms of dust attenuation by Bu at et ai] (l2007a) 
(see also lBurgarella et alj d2007l) ). 

Some galaxies of the IR selection, detected in UV, do not 
appear in the UV-selected sample. The same is true for some 
galaxies of the UV selection that have an IR detection but are 




10 11 

log(Lb i) (solar units) 



Fig. 1. IR to UV luminosity ratio as a function of the bolometric 
luminosity for both samples at z=0 and z=0.7. At z=0.7 individ- 
ual data are plotted: UV selection with blue "plus" for galaxies 
also detected at 24//m and arrows down for galaxies not detected 
at 24/zm; IR selection with red crosses for galaxies also detected 
at UV and arrows up f or galaxies not de tected at UV. At z=0 
volume average trends dBuat et al.Ll2007bl) are plotted with solid 
lines and error bars: the UV selection with a light solid line and 
the IR selection with a heavy solid line 



not included in the IR-selected sample: these "single" objects 
appear in Fig.[TJas simple crosses (IR selection) or "plus" sym- 
bols (UV selection) and represent 30% of both samples. Their 
presence is due to the method used to cross-correlate the data: 
for the IR (resp. UV) selection, the cross-correlation with the 
UV (resp. IR) sources is made using the complete list of the UV 
(resp. IR) detections before their identification with COMBO 17 
objects and any other selection. The truncations to have com- 
plete volume-limited samples account for a half of the "single" 
objects. The other half of the IR (resp UV) "single" objects have 
a UV (resp IR) counterpart that has not been identified with a 
COMBO 17 source. Such a rat e of non-identification ( 15 %) is 
consistent with those quoted by Burga reUa et alj d2007l) . 

3.2. Color distributions 

It is well established that a population selected in optical or in 
NIR exhibits a bimodal distribution of color s. This result, first 
obtained at low z (e.g . iBaldrv et all 120041) . is also observed 
at higher z (e.g iBellj. 12004 ICooperet all 120071: iBell et all 
120071 : lElbaz et alll2007l) . Blue galaxies are star-forming objects, 
whereas the redder ones are quiescent systems. Our galaxies 
were selected either in UV or in IR and so are likely to show 
active star formation. For the sake of comparison with other 
studies, we can check what sort of galaxies we are dealing with 
in terms of optical colors. Unfortunately for us, this check is 
very difficult to perform at z=0 since our sample was built from 
2000 deg 2 not covered by homogeneous optical surveys like 
the SDSS. Nevertheless, we can rely on other studies of UV or 
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IR selected galaxies. Iglesias-Par amo et alj d2007l) studied UV- 
selected galaxies from z=0.2 to z=0.7. The U-V distribution of 
their sample remains unimodal over the full redshift range they 
analyzed: according to the classical subdivision into blue and 
red galaxies , only blue star-forming galaxies are selected in UV. 
iGotol d2005l) studied the optical properties of a sample of galax- 
ies selected from the IRAS catalog. The color distribution that 
he found is broad and extends from the blue peak to the red one 
defined from the SDSS studies. 

At z=0.7 the situation is more favorable to our galaxy sam- 
ples. We can work with the R-I color from the EIS catalog, which 
is similar to the u-g color in the rest-frame of the galaxies. The 
R-I color is conve rted in AB magnitude s according to the conver- 
sion formulae of lArnouts et alj d2001). The color distributions 
are displayed in Fig. [2] The distinction be tween red and b lue 
populations is below and above R - I =2 1 dElbaz et aUl2007l) . 

The red population is under-represented for both selections 
but the two distributions clearly differ. As expected in UV we se- 
lect preferentially blue galaxies that are active in st ar-formation. 
The se lection looks similar to what is found by Bun dv et alj 
(2006) for the DEEP2 survey restricted to galaxies with a SFR 
larger than 0.2 M yr _1 measured with the [Oil] equivalent 
width. The distribution of the IR-selected galaxies is broader 
and shows a tail towards redder objects, but we do not see any 
clear bimodality as reported for example for optically selected 
galaxies (see refe rences above). This agrees with the findings of 
iBell et ail d2005l) at a similar redshift that the galaxies detected 
at 24 /urn are spread over a wide range of o ptical colors . If we 
compare, for example, to the distribution of Elbaz et al. (2007) 
it seems that our distribution is broader with a substantial frac- 
tion of our sample between the two peaks defined in U-B rest- 
frame (similar to observed R-I at z=0.7). Therefore a selection 
at 24 jt/m does not exhibit bimodality, and it also seems to be the 
case at z=0 for galaxies selected from IRAS: IR-selected galax- 
ies seem to populate the "green valley" located between the blue 
and red peaks. A substantial fraction of these galaxies with in- 
termediate colors exhibit a strong dust attenuation with the IR to 
UV flux ratio greater than 10 and they are often not detected in 
UV (cf Fig. [1). It may be noted that the fraction of non detec- 
tions does not vary a lot with the bolometric luminosity of the 
galaxies. The IR selection does not include the large number of 
very blue galaxies detected in the UV with low bolometric lu- 
minosity (cf Fig. Q]). Last of all, when only galaxies selected at 
both wavelengths are concerned (dashed histograms) the prop- 
erties of both samples are found to be similar, it is also clearly 
seen in Fig.Q] 



4. Specific star formation rates 

4.1. Variation in the SSFR within each sample 

In this section, we analyze the variations in the specific star for- 
mation rates (SSFR, defined as the total SFR divided by the stel- 
lar mass) as a function of the stellar mass within our samples of 
star-forming galaxies. In Fig.[3]the SSFRs versus stellar masses 
are plotted for both samples at bot h redshifts . Averag e trends are 
plotted at z=0. They come from iBuat et all d2007bl) where vol- 
ume averages were performed. The values were converted using 
the SFR and stellar mass calibrations adopted in this paper (sec- 
tion 2.1). At z=0.7, the values for each galaxy are plotted. 

Before any interpretation, we must estimate the detection 
limits in SSFR to determine where the results are reliable, on the 
basis of the limiting fluxes adopted in UV and IR. Stellar masses 
were estimated from the 3.6 yum flux, available for each galaxy, 
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Fig. 2. Observed R-I color (AB scale). Heavy solid line: whole 
IR selection, heavy dashed line: IR selected galaxies also de- 
tected in UV. Light solid line; whole UV selection, light dashed 
line: UV selected galaxies also detected at 24 jum 



and are thus unaffected by detection limits. Practically speaking, 
for the IR selection we adopt the limit of log(LiR) = 10.34 (L ) 
used to build the sample (cf section 2.2.1) and translated it into 
SFR at z=0.7 with Eq. 1. We obtain SFR hm = 2.34 Movr -1 . 
Given the high values of the IR to UV ratio for the IR selection 
(cf Fig. [TJ, we have neglected the contribution of the UV emis- 
sion to estimating SFR\i m . For the UV selection, the adopted 
limit log(Luv) = 9.39(L ) translates into an "unobscured" SFR 
at z=0.7 of 0.50 M Q yr _1 according to Eq. 2. However we must 
also account for the contribution of the IR emission to the SFR 
which is not negligible at all (cf Fig. [TJ. We used the detection 
limit of 0.016 mJy at 24 yum (cf. section 2.2.1) translated into 
an "obscured" SFR at z=0.7 of 1.15 M Q yr _I . The limit thus ob- 
tained for the total SFR is SFR\ lm = 1.65M yr~'. The resulting 
detection limits obtained for z=0.7 for the SSFRs (and only con- 
strained by those on the SFR) are reported in Fig. [3] 

The z=0 samples are purely flux-limited so we cannot es- 
timate t he detectio n limits in the same way as for the z=0.7 
sample. IBuat et all (l2007bl) built bolometric luminosity func- 
tions with these samples (Lboi = ^uv + (1 — if) ^ir) down to 
6.3 10 s L and 1.8 10 9 L for the UV and IR selections re- 
spectively. These luminosities were used to estimate the limits 
reached in SFR and they can be translated in SFR using Eq. 
1. We obtain SFR > O.O7M yr _1 for the UV selection and 
SFR > O.^Meyr 1 for the IR selection. Such low SFRs put 
detection limits well below the values reported in Fig. [3]for the 
mean trends at z=0. 

The influence of the choice of 77 (Eq. 3) can be easily 
checked: the measure of the SFR is mostly dominated by the 
IR emission with only a small contribution from the UV one, es- 
pecially for massive galaxies which will be the main topic below. 
Therefore adopting 77 = 0.3 instead of r\ — at z=0.7 will reduce 
the SFR and hence the SSFR by a factor 0.7, i.e. a translation of 
-0.15 dex on the vertical axis of Fig. [3] 
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Fig. 3. SSFR versus stellar masses. Left panel: IR selection, the z=0.7 galaxies detected at 24 fim and in UV are plotted with red 
crosses, the galaxies only detected at 24 fj.m with dots. The solid line with error bars (1 cr) is the result of the volume average at 
z=0.7, the dashed line the volume average at z=0. The diagonal line represents the detection limit at z=0.7. Right panel: Same plot 
for the UV selection. At z=0.7 galaxies detected at 24 yum and in UV are plotted with blue crosses, the galaxies only detected at UV 
with black crosses, the lines are defined in the same way as for the IR selection. Adopting 77=0 instead of 77=0 at z=0.7 would shift 
all the data at z=0.7 by -0. 15 dex along the vertical axis 



It is obvious from Fig. [3] that the locus of the galaxies with 
masses lower than 10 10 M is completely governed by detection 
limits at z=0.7. The observed trends are only reliable for M > 
10 10 Mq. In the following, we limit our analysis to this mass 
range. The trends found for both selections seem to be similar, 
they will be compared in the next sections. 

Our samples were built to be complete in terms of the star 
formation rate (as discussed above in this section) but are not 
expected to be complete in mass. Indeed a mass-selected sam- 
ple also contains quiescent galaxies that are not present in our 
sample. At z=0 we built the stellar mass function using the 
V IV max formalism for the IR and the UV-selected sam ples. They 
are co mpared to the stellar mass functions derived by Bel l et al] 
d2003l) from SDSS +2MASS data. W e have applied the correc- 
tion preconised by iBell et ail d2003l) for a Kroupa IMF. As ex- 
pected, the stellar mass distributi on of the UV-s elected sample 
is similar to what is found by IBell et ail d2003l) for late-type 
galaxies (selected according to their colors), whereas the IR se- 
lection leads to a larger number of massive galaxies but still 
below the total stellar mass distribution. At z=0.7 we can di- 
rectly compare the counts obtained for a selection at 3.6 microns 
(IRAC), 24 microns (MIPS), and in NUV (GALEX) for the red- 
shift bin 0.6-0.8. Fig. [5] summarizes the comparison where the 
counts are reported against the stellar mass. The selection at 
3.6 microns consists of sources brighter than 3 /vJy, which cor- 
responds to 10 10 M Q at z=0.7; this limit is very conservative 
given the depth of the IRAC observations dElbaz et al.L l2007t 
ISanders et all 120071) . Both MIPS and GALEX selections lead 
to lower counts than the IRAC one. The effect is more extreme 



for the UV selection. This was expected from Fig Q] where it is 
obvious that the UV selection misses bright (and therefore mas- 
sive) galaxies with a substantial dust attenuation. Conversely the 
contribution of the UV sources becomes significant for masses 
around 10 M . We have also included in Fig. [5] galaxies se- 
lected in either the UV or the IR in order to analyze their total 
mass distribution. The difference between the number of galax- 
ies in our GALEX+MIPS selections, and the IRAC one is likely 
to be due to the presence of quiescent systems detected neither 
at 24 yum nor in NUV. It is qualitatively consistent with the rela- 
tive contribution of spira ls and irregu l ars to the total stellar mass 
function as measured by Bun dv etaTI (|2005) in the same redshift 
range. A full analyzis of the luminosity functions in our selected 
samples is in preparation (Takeuchi et al.). 

Before comparing to models we also calculate average trends 
at z=0.7 (Fig. |3J. There is no need for volume corrections (cf 
section 2.2) and simple mean values can be calculated. We have 
used the Kaplan-Meier estimator for the IR selection and the er- 
ror bars plotted in Fig.[3]correspond to the 1 cr dispersion. For the 
UV selection, the presence of upper limits located on one side of 
the distribution hampers the use of the Kaplan-Meier estimator. 
Therefore we have calculated the median instead of the mean for 
the UV selection at z=0.7. Very similar trends are found in both 
selections (see also Fig. [7] where average trends for both selec- 
tions are overplotted). Even if the UV selection misses highly 
reddened galaxies (cf Fig. [TJ, the star formation activity of the 
galaxies selected in UV is found to be similar to those selected in 
IR: the SSFRs do not depend critically on the adopted selection 
(as long as they are estimated by combining UV and IR emis- 



V. Buat et al.: Star formation history of galaxies from z=0 to z=0.7 



7 



T3 

ra 



M 
O 



- 

- =r ~~ — 

£ -~ ■ 




1 | 1 1 1 1 | 1 1 


_ 

_ 

— 


^ o-^l"^ — ~~ 






: 








- 










- 






- 






\ \ \ \ 




1 




X \\ \ 


~- 






V /M 








\ \/ \ 








Y\ \ 








V N \ 
\ \ 








\ 




" 










8.5 9 9.5 10 10.5 11 11.5 
lo g( M star) (solar units) 



Fig. 4. Stellar mass functions for the z=0 samples selected in IR 
(circles) and in UV (triangles). The solid l ine withou t symbols is 
the total stellar mass function derived by iBell et all d2003l) and 
the dashed one denotes the stellar mass function for late type 
galaxies by the same authors 
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Fig. 5. Galaxy counts as a function of the stellar mass for a se- 
lection at 3.6 /im (stars), 24 /im (circles), in NUV (triangles) and 
for the galaxies selected either at 24 fim or in NUV (squares) 



sions), a result already found a t z-0 once volum e corrections 
are applied to the local samples dBuat et aUl2007bl and Fig.[7]>. 



4.2. Comparison with previous studies at intermediate 
redshift 

Numerous works address the question of the variation in the spe- 
cific star formation rate at various z. Whereas these studies con- 
clude to a general decrease of the SSFR with galaxy mass even at 
high z (the so-called downsizing effect), the trends exhibit large 
departures from one work to another. Our aim in this section is 
not to perform an exhaustive comparison with all the available 
results: we will instead discuss a few studies performed in the 
same redshift range as the present work, and using methods to 
estim ate SFR and stellar masses similar to ours. 

Zhen g et all (|2007a) analyzed the dependence of star for- 
mation on galaxy mass within the COMBO 17 survey adding 
SPITZER data to measure accurate SFRs. The SSFRs they re- 
ported are consistent with ours within error bars but with a 
steeper decrease in the SSFR when stellar mass increases com- 
pare d to what we find . Their sample contains blue and red galax- 
ies, iBeil et alj d2007l) considered them separately and found a 
flatter distribution of SSFR for the blue ones than for the red 
ones. In a similar wav lElbaz et al.l d2007l) analyzed the GOODS 
survey and found a flat distribution of SSFRs at z=0. 8-1.2 and 
between 10 10 and 10 n M o . The decrease in the SSFR they found 
around lO n M is clearly due to the red galaxies, whereas blue 
gala xies exhibit a rough ly constant SSFR (their figure 17). 

iNoeske et al.1 d2007l) also found a flat distribution that is 
fully consistent with ours within the AEGIS survey (both in 
absolute values and trends) for galaxies with stellar masses be- 
tween 10 10 and 1O U M and more massive galaxies have a lower 
SSFR. Studies based on GALEX data and combining UV and 
IR emissi ons have also led t o an almost flat distribution of 
SSFRs. Iglesia s-Paramo et al.l d2007l) performed a UV selection 
combined with SWIRE data from z=0.2 to z=0.7 and found 
a flatt ening of the SSFR v ersus M stal variation from z=0.2 to 
z=0.7. iMartin et all d2007l) also found a flat distribution in the 
CDFS at intermediat e z with deep GALEX and MIPS data. 
Zamoiski et al.l d2007l) perform a very complete analysis of the 
COSMOS field (although only based on UV-optical rest-frame 
data) at z=0.7. Their UV-detected sample also exhibits a rather 
flat distribution of SSFR, much flatter than obtained for their en- 
tire sample (including objects not seen in UV). 

Although a large dispersion is found in the already pub- 
lished results, we find good agreement with those focusing on 
star-forming galaxies alone. The steeper decrease in the SSFR 
with increasing mass found in some studies at intermediate 
redshift is likely to be due to the presence of quiesce nt systems, 
which are absen t in o ur present selection (e.g. IBell et all 
120071: lElbaz et al.L 120071) . Similar differences are found at z=0 
when only star- forming galaxies are se l ected or quiescent 
are also added dBrinchmann et all l2004t iBuat et all l2007bt 
Elba z et ail 120071) . Most of the decrease found at z=0.7 also 
occurs for masses hig her than 10 u M n (eg. lZheng et al.l [2007ah 
INoeske et all 120071: lElbaz et all 120071) which are not well 
represented in our sample. 



5. The models 

Recent studies have found evidence of the minor role of strong 
mergers in the evolution of massive galaxies from z=0 to ~ 1 
dBell et all 12001 iMelbourne et all 12051: IZheng et all l2007bl) . 
Therefore it is tempting to try to fit the data with models that 
assume a smooth average evolution, which is weakly challenged 
by small interactions and minor mergers. 
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Table 2. Models with A = 0.05, Zf - 6, stellar masses are calculated at z=0 (column 2) and z=0.7 (column 3) with the analytical fit 



velocity (km s 1 ) 


log(M star ) c =o(M )) 


log(M sl!U .) c =o. 7 (M ) 


a 


b 


c 


80 


8.89 


8.19 


6.62 


0.41 


0.36 


150 


9.92 


9.51 


8.74 


0.98 


-0.20 


220 


10.52 


10.25 


10.01 


1.25 


-0.55 


290 


10.94 


10.75 


10.81 


1.35 


-0.74 


360 


11.25 


11.10 


11.35 


1.37 


-0.85 



Recently Noe ske et al] (|2007) proposed a simple model of 
gas exhaustion to interpret the variation in the SFR as a function 
of the stellar mass out to z=l.l. In their model galaxies exper- 
iment exponential star formation histories with e-folding rates 
and redshift formation varying with the galaxy mass, less mas- 
sive galaxies being younger and having higher e-folding rates. 

Here, instead of building ad hoc models (the parameters in 
INoeske et all [2007, are fine-tuned for t heir data) for our study , 
we follow a "backward" approach (e.g. lSilk & Bouwensl [l999) 
by starting from studies of the Milky way and the nearby uni- 
verse to extrapolate the behavior of galaxies at higher z, without 
any further adjustment of the models that are constrained in the 
local universe (see details below). 

5.1. Description of the models 

To interpret our results with physically motivated models, we 
use a grid of models for the evolution of spiral galaxies that 
is similar to the one presented in Boissier & Prantzos (2000). 
These models were calibrated to reproduce many properties of 
the Milky Way (Boissier & Prantzos, 1999), and were compared 
to nearby spirals in subsequent works. They should be adequate 
to broadly represent the family of spirals from very massive ones 
to irregular low-mass exponential disks. The models use scaling 
relationships to simulate disks of various rotational velocity and 
spin para meter (A, measuring the specific angular momentum, 
see e.g. iMo. Mao & Whitel [T998, for a definition). The veloc- 
ity is closely related to the total mass of the galaxy and, as a 
result, its final stellar mass. Each stellar mass (at any redshift) 
corresponds to one and only one model for a given specific an- 
gular momentum A. Th e specific angular momen tum has a log 
normal distribution (e.g. lMo. Mao & WhiteLll998l) . Spirals have 
spin parameters between 0.02 and 0.08, with a typical value 
of 0.05 (the Milky Way G alaxy corresponding to A ~ 0.03 in 
iBoissier & Prantzoa, l2000h . Since we are interested in general 
trends with mass, we only use models with this average A =0.05 
in the following. The quantities we present in this paper present, 
in any case, a much stronger dependence on the velocity than on 
the spin parameter. 

The models assume that galaxies are formed by progres- 
sive infall of primordial gas, starting at high redshift (z/ = 6). 
Models of the Milky Way shown early-on that infall of "fresh" 
gas provides a good explanation especially for the distribution 
of metallicity of G-Dwarfs in the Solar neighbo r hood and the 
idea has been widely used (e.g. Larsonl fl972k iPagell 119971; 
Chiappini et allfl997HBoissier & Prantzoslll999l) . Such models 
do not state that this gas has been forever in a reservoir around 
the galaxy. In a more modern context, it could very well be that 
this infall corresponds to small satellites (with large gas frac- 
tion) being accreted by the galaxy during minor mergers. Inside 
the disk, stars form from the gas following a Sc hmidt-like law , 
includ ing a dynamical factor. In a recent work (Bo issier et all 
120031) . the star formation law was empirically determined to be 

(4) 



relating the surface density of star formation rate £sfr 
(M pc _2 Gyr _1 ) to the gas surface densities Zgas (M pc~ 2 ). Here 
V(R) is the rotation velocity (km s _1 ) at radius R (kpc). This for- 
mulation is very close to the one originally used, and subsequent 
models (includi ng in this paper) have used it . The results are very 
close to those in Boissier & Prantzos (2000), and show the same 
global trends. 

The star formation history in these models strongly depends 
on the built-in assumption that infall proceeds at a higher pace 
early in the history of massive galaxies with respect to the Milky 
Way, and later in lower mass galaxies. This assumption corre- 
sponds to the so-called "downsizing" effect, and was made to 
adequately reproduce observed trends in nearby spirals. 

By construction, these models should match the z=0 uni- 
verse at the current epoch. The use of such models in a backward 
approach is indeed justified only because a very large number of 
properties of nearby galaxies (i.e. the redshift zero universe) are 
correctly reproduced by them. It is notably the case (and we re- 
fer to the references for further details) of scaling relationships 
such as surface brightness and scalelength versus magnitude, the 
Tully-Fisher relationship (and its dependence on wavelength), 
color-magnitude diagr ams (B-K vs K), lumino sity-metallicity 
relationship, spectra dBoissier & Prantzoa |2000|). colors and 
abundance gradients in spirals (Prantzos & Boissier, 2000), star 
formation rates and gas fractions (Boi ssier et all 120011) . 

In a "backward" approach, their prediction about the past 
was only compared to a few properties of small samples of spi- 
rals at higher redshift dBoissier & Prantzosl 1200 lb . and no seri- 
ous discrepancies were found with the data available at the time. 
Note that the observed B-band luminosity function at redshift 
zero was reproduced by the models in this paper. This result was 
obtained by construction since the authors assumed a circular ve- 
locity distribution Vc derived from observations, and the models 
reproduced the Mb vs Vc relationship. 

We propose here to use these models to interpret data con- 
cerning the SSFR (such data did not exist a few years ago). 

5.2. Star formation histories 

The star formation histories (presented in Fig. [6] left panel) of 
the models result from the non linear combination of a star for- 
mation law, an infall history and its mass dependence. To provide 
an easy way to compare our star formation histories to others, or 
to various data sets, we propose an analytical form for the star 
formation history of our models. After testing several possibili- 
ties, we obtained a correct fit using three parameters: 



SFR(t) = Af fc exp(-(f/r) as ) or 
log(SFR(t)) = a + b log(f) + c f 5 
with a = log A c= -0.43r~ a5 . 



(5) 
(6) 
(7) 



-'SFR 



2.63 10- 3 2^ 8 S V (R)/R 



We can assign a tentative physical significance to each of 
the parameters : A is a scaling factor, b is measuring the rate at 
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Fig. 6. Left: Star formation histories of the adopted models for five final stellar masses (equivalently: five rotational velocities), from 
the lowest mass (bottom) to the highest one (top). The dotted line is the result of computing of the chemical evolution models and 
the solid line is the analytical fit of (Eq. 6). On the x-axis, T=0 Gyr corresponds to the formation of the galaxies (Z/ = 6), z=0.7 
corresponds to T = 6.2 Gyr and z=0 to T = 12.5 Gyr. Right: Variation in the three parameters (a,b,c) of the analytical star formation 
history (Eq. 6) as a function of the present (z=0) stellar mass (values are given in Table 2 for the five models shown on the left, also 
shown as dots in this panel). 



which the star formation increases early-on, and t is a time-scale 
for the decay of star formation at later time (this is not always the 
case: for low-mass galaxies the star formation rate is still rising 
at late epochs). The values of these parameters for the A = 0.05 
models are given in Table 2 for a few values of the circular ve- 
locity and their dependence on the final (z=0) stellar mass is 
shown in Fig.[6](right panel). The quality of the fit can be judged 
from the left panel. With such a formula, it is straightforward 
to compute analytically the evolution of stellar masses assum- 
ing the instantaneous recycling approximation (e.g. P agell[T997[) 
and a returned fraction of R=0. 3 for the iRroupa et al.l (1993) 
IMF dBoissier & Prantzosl 1 1999b . We checked that this simple 
approach gives a very good approximation to the results of the 
numerical computations implementing infall, star formation law, 
and finite lifetimes of stars. 

Before comparing the models and the data, we should note 
that some ingredients of the chemical evolution models are 
rather uncertain. Especially, the star formation law efficiency 
is very dispersed among disk galaxies, making the star forma- 
tion rate uncertain b y about a factor 2. In addition, we used the 
iKroupa et al.l (l993) IMF, for consistency with earli e r work s. We 
checked that if we used the more recent iKroupal ([2001) IMF 
for the models, the star formation rate and stellar mass during 
the overall history of the galax y would change by on ly 10 to 
20 % Modifying the IMF from lKroupa et alj J 19931) to lKroupal 
(2001) mainly affects the population of massive stars and, as a 
consequence, the calibration of the current star formation rate 
and the predict ed UV fluxes is why we use the more recent IMF 
dKroupall200ll) to derive SFR from the UV fluxes (Eq. 2). Model 
predictions other than UV fluxes depend much more weakly on 
this choice. 

iHammer et al.l ([2007) suggest that the Milky Way, hav- 
ing escaped significant merging over the last ~ 10 Gyr, is 
in fact unrepresentative of spirals in contrary to M31. Our 
models, even if they were calibrated in the Milky Way, suc- 
cessfully reproduced many properties o f nearby spirals, in- 
cluding abundance and colour gradients dBoissier & Prantzosl, 



l200Ct iPrantzos & Boissierl l2000h . The star formation law (the 
"Schmidf'-like law mentio ned above) was dete rmined from 16 
galaxies with Ha profiles (Boiss ier et all 12003). and is consis- 
tent with the SFR indicators in the Milky Way. The SFR pro- 
files have been deter mined in 43 nearby spirals from UV data 
dBoissier et al.L 120061) and compared to gas profiles. The SFR- 
gas connection they find agrees with the one derived from Ha 
profiles, even if the relationship extends to lower gas surface 
densities in the UV (the untypical galaxy with this respect was 
actually M31, which does not follow the Schmidt law in its in- 
ner part). We recognize that the models calibrated in a relatively 
calm Milky Way probably not include all the details of disk 
galaxy evolution. However, they do agree with many properties 
in nearby spirals as discussed above and we believe they repro- 
duce the major traits of disk evolution even if real galaxies are 
likely to sometimes suffer interactions not taken into account in 
this approach (as long as they are not dramatic events destroying 
the disks such as major mergers. 

6. Results and discussion 

6. 1 . Comparison between models and data at z=0 and 
z=0.7 

In Fig. [7] the results of the model are compared to the mean 
observational trends obtained in section 4. The agreement can be 
considered as good for both selections at M > 10 10 M given the 
large intrinsic dispersion of the data and the fact that no fitting 
has been performed. At z=0.7 the agreement is remarkably 
good given the uncertainties inherent to both the models and the 
data (see the discussion below). At z=0 the SSFRs predicted by 
the model lie slightly above most of the mean observed values 
and exhibit a steeper decrease as the stellar mass increases, as 
compared to the UV selection, although the intrinsic dispersion 
of the observed quantities is very large. As discussed in section 

5.2, models themselves contain some uncertainties. Modifying 
the star formation law efficiency within the uncertainties could 
modify the star formation rate by a factor ~2, and various IMF 
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may modify the results of the models by up to ~ 20 % as seen 
above. The models errors are likely to be systematic if the IMF 
and the star formation are universal, i.e. directly linked to the 
physics of star formation on a local scale rather than depending 
on the global properties of galaxies. In that case, the trends 
obtained by the models are robust even if the absolute values 
could change. 

Models predict directly the stellar masses, their estimates do 
not depend on any assumed mass-to-light ratio. The uncertainty 
on the mass-to-light ratios only affects t he masses derived from 
the ob servations. It can reach a factor 2 dKannappan & GawiserL 
120071) . similar to the dispersion found for the observed SSFRs 
for a given range of stellar mass. Therefore both models and 
observables are likely to be affected by systematic errors large 
enough to explain the slight shift found at z=0 between the mod- 
els and the data in Fig. [7] Another issue may be that the SSFR 
trend with the stellar mass is flatter for the UV-selected sample 
than the models predictions. Indeed the UV selection leads to 
a selection of late-type systems (cf section 4.1 and Fig.|4]i; and 
since the contribution of early type galaxies increases with the 
stellar mass (Fig.|4|, we expect a flatter distribution of the SSFR 
for the late type systems alone as compared to the total galaxy 
population. The models are set up to reproduce the mean evolu- 
tion of galactic disks. For massive galaxies, they correspond to 
relatively early disk types (as early as Sa) in which the specific 
star formation rate is low: i.e., the star formation activity is low. 

The UV selection can pick up massive galaxies with ac- 
tive star formation due to burst or refueling through interactions 
(leading to later type galaxies for the same mass) not introduced 
in the models that then correspond better to the galaxies from 
the IR selection (picking up massive galaxies with more quies- 
cent star formation). 

At higher z, because of the higher star formation activity in 
all the galactic disks, this effect is expected to be smaller and the 
dis crepancy between the models and the data disappears. 

Noes keet al.l ((2007) have performed a similar analysis to 
ours. They selected galaxies between z=0 to z=l from the 
AEGIS survey. They derived SFR by combining IR, optical, and 
UV data. Their SSFR distributions are very consistent with ours. 
As a result, the gas exhaustion model they propose would also 
fit our data reasonably well. In their model, star formation rates 
follow an exponential decay after a formation redshift Zf. Less 
massive galaxies have a longer e-folding rate and a lower redshift 
of formation. Nevertheless, it implies a wide range of redshift 
formation, from z/=l for 10 K) M G up to Z/=3 or lO n M . In our 
model, on the contrary, all the galaxies have the same formation 
redshift zj = 6, which is the epoch when galactic building blocks 
are assumed to start to exist. It is rough ly the redshift w hen the 
first galaxies are confirmed to exist (e.g. lSchaerert 120071) . In this 
approach, it is the time variation of the SFR that depends on the 
galaxy mass, in such a way that the bulk of star formation oc- 
curs at different ages according to the galaxy mass. Actually, the 
star formation histories obtained in our models a r e qu alitatively 
similar to the schematics proposed by Sandag e dl986l). which 
inspired more recent works. Among them, Gavazzi etjij] d2002l) 
used a sta r formation histo ry law "a la Sandage", which is mim- 
icking the Sandage ( 1986) trends. Quantitatively, our star forma- 
tion history is slightly different: for massive galaxies especially, 
the SFR does not decrease as quickly at late epochs. 

Because of this, the mean evolution of the SSFR from z=0 
to z=0.7 in star-forming galaxies of 10 10 - 10 u M o is consis- 
tent with the one predicted by a simple but physically motivated 
model of secular evolution, in which the galaxies are progres- 
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Fig. 7. Specific SFRs versus the stellar mass: comparison be- 
tween data and model. UV selection: blue circles and error. IR 
selection:red triangles and error bars. The filled symbols are for 
z=0, the empty ones for z=0.7. The solid lines represent the 
model summarized in Table 1 (see text) 



sively built by accretion of low-metallicity gas and in which star 
forms according to a "universal" law of star formation. It does 
not mean that individual galaxies cannot experience any burst 
events or interactions, but our aim is only to reproduce mean 
trends. At z=0, act ive galaxies were found in the IR selection 
dBuat et all l2007bl) . but these objects were also diluted among 
the more numerous fainter objects when volume averages were 
performed. In the same way the most active galaxies of our sam- 
ples at z=0.7 have SFRs reaching 40 M yr , a factor ~ 2 larger 
than highest SFR from our models (Fig. [6]). Indeed luminous in- 
frared galaxies (LIRGs which account fo r 30 % of our IR selec - 
tion) are known to experience bursts (e.g. Marcillac et all l2006). 
Nevertheless when averages are calculated their influence is di- 
luted. Indeed, our results do not imply a secular evolution for the 
whole galaxy population but are instead consistent with a smooth 
mean evolution in agreement with the idea that galaxies are built 
by infall (which can be in the form of minor mergers, as sug- 
gested above). In this context, major mergers or other extreme 
events are not expected to participate much in the construction 
of the average galaxy. 

Interestingly, the model with a circular velocity of 220 km/s, 
corresponding to the one of the Milky Way, has a stellar mass of 
log(M star ) = 10.52 (Mo) at z=0 and log(M star ) = 10.25 (Mg) at 
z=0.7. In other words, our samples probe Milky Way analogues 
in terms of stellar mass. 



6.2. Comparison with semi-analytical models 

The models presented above have the advantage of reproduc- 
ing many properties of galaxies at redshift zero. They assume 
smooth evolutions on average, which might not be realistic in 
the paradigm of hierarchical galaxy formation. Semi-analytical 
models (SAMs) usually follow a similar approach to our mod- 
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els (e.g. assuming empirical recipes for star formation laws and 
scaling relationships) but obtain the mass accretion and merging 
histories of the galaxies by following the hierarchical growth of 
dark haloes (and the baryonic galaxies within them). 

We also compared our data to the results of two sets of 
SAMs. iKitzbichler & White! d2007l) have produced lightcones 
simulations that can be very convenient ly compared to ob- 
servations. They applied the SAM of Croton et al. (2006) 
to the Millennium R un simulation ( Spring e! et all 120051) . 
Nagashim a et al.l ((2005) have constructed a numerical catalog 
also based on SAMs and combined with high resolution N-body 
simulations. We have applied selections to these simulated cata- 
logs as close as possible to those obtained for our observational 
data (section 4.1): an SFR higher than 2 M yr~' for galaxies 
at a redshift between z=0.6 and 0.8 and an SFR higher than 0.1 
M© yr~' for nearby galaxies (z < 0.1). We computed average 
trends with the simulated data in the same way as in section 
4 for our observational datasets. The results are displayed in 
Fig. [8j where the simulated quantities are compared to our 
mean observational trends obtained for the IR and UV selected 
samples. It can be seen that the simulated data also reproduce 
the values of the mean SSFR observed at z=0.7, although the 
simulations lead to a somewhat steeper decrease in the mean 
SSFR towards high stellar masses than observed. A rather large 
discrepancy is observed at z=0 for the millennium simulations: 
the simulated galaxies have lower SSFRs than the values 
measured in our samples. If we ignore this discrepancy the 
agreement is satisfactory. We can conclude that SAMs present a 
similar mean behavior to the one found for our simpler evolution 
model aimed at reproducing present day spirals without strong 
merger events. This result reconciles backward and forward 
approaches for the evolution of galaxies at least until z=0.7 and 
validates our simple approach for star-forming galaxies. 

The SAMs should be able to reproduce galaxies with very 
high SSFR induced by major mergers, whereas our model only 
aims at reproducing mean trends. Nevertheless, it may be mean- 
ingful to mention that the timestep used in the simulation (about 
300 M yr for the Millennium simulation (IKitzbichler & White, 
120071) ) might be too large to catch short time bursting systems 
and the SFRs averaged over such a time step cannot be very 
high. This aspect of the simulation works should be improved 
in order to make a more extensive comparison of star-formation 
related quantities, like SSFR. A more exhaustive discussion of 
the results of semi-analytical simulations is beyond the scope of 
this paper and will be shown elsewhere. 



6.3. Predictions and comparison with data for higher z 

Since our models are able to reproduce the evolution of massive 
star-forming galaxies from z=0 to z=0.7, we can take a step fur- 
ther and make some predictions at higher z. These predictions 
can be considered as a reference for quantifying the expected 
properties of distant galaxies in the absence of major merging, 
implying a strong starburst and a very high efficiency of star for- 
mation. In Fig. [9] are reported the evolution up to z=4 of the 
SSFR and of the gaseous phase metallicity of the chemical evo- 
lution models. 



6.3.1 . Redshift evolution of the specific star formation rate 

A main prediction issuing from such a secular evolution is the 
reduction in the range of values of the SSFR when the redshift 
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Fig. 8. Predictions of the Millennium model at z=0.7 and 
z-0 (dot-dashed lines ) and of the numerical catalog of 
Nagashim a et al.l d2005) (dashed lines). All the other lines and 
symbols are the same as in Fig [7] The mean error bar (lcr) is 
indicated on the right side of the panel 



increases (Fig.|9j panel a). This is a natural result because for a 
hypothetical first generation of stars formed during a brief time 
6t, we should have S FR x 5t ~ M stm -. Equivalently the model 
predicts a flattening in the variation of the SSFR as a function of 
the stellar mass as z increases (Fig. [9] panel b). Such a flattening 
is indeed observed in our data up to z=l and i n other studies 
(llglesias-Paramo et all [20061: iMartin et al.Ll2007l) 

At high z, very few observational studies exist. Moreover 
they are often difficult to compare because selection effects are 
likely to become extremely strong. For e xample a K se l ection 
at z=4 corresponds to a B rest-frame one. Feulneretal. (2005) 
analyzed the FORS Deep Field and GOODS-S field to derive 
the variation in the SSFR up to z=4 for different ranges of stel- 
lar mass from an I and a K selection. Their average values for 
intermediate mass systems (10 9 5 - 10 10 5 M o ) are in reasonable 
agreement with ou r models from z=l t o z=3. In all the redshift 
range they explore. IFeulner et al.l (|2005) find a strong decrease in 
the SSFR as the mass increases, which might seem at odds with 
the predictions of o ur model (Fig. [9] panels a & b). However, 
IFeulner et al.l d2005l) did not select only star-forming galaxies, 
and the contribution of quiescent objects can steepen the varia- 
tion in the SSFR with the mass since they are essentially very 
massive objects. Indeed very recently Da ddi et al.l (|2007) stud- 
ied star formation in massive galaxies forming stars actively up 
to z~2 and found a roughly constant SSFR at z=2 for galaxies 
detected at 24 /mi and with a mass comprised between ~ 5 10 9 
and ~ 10 11 M , in agreement with the predictions of our model. 
iDaddi et alJ 0007) find an increase in the SFR at a given mass 
from z=0 to z=l and z=2 by a factor ~4 between z=l and 2 and 
a factor ~30 between z=0 and 2, these values are higher than 
those we predict (a factor ~2.5 between z=l and 2 and a factor 
~10 between z=0 and 2). 
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lLamareille et al.l d2007l) present the variation in the SSFR up 
to z=2.5 for the VIMOS VLT Deep Survey in different bins of 
stellar mass. The absolute values, as well as the variation they 
found for the galaxies with 10 < log(M star ) < 11, are very con- 
sistent with our models. From their Fig 2 there is also a hint of 
some flattening of the SSFR-M star variation when z increases but 
still of a decrease i n the SSFR as the m ass increases. As in the 
case of the study of lFeulner et alj d2005l) . the VVDS is an I band 
selected survey that might not only select star-forming galaxies 
but also more quiescent objects: a steep decrease in the SSFR 
with increasing stellar mass is expected with a bimodal galaxy 
population. 

SAMs have been used in the literature to compare to 
observation a l sets in the same way we proceeded at z=0.7. 
iDaddi et all d2007l) simulated star-forming galaxies at z=2 from 
the Millennium database. The regression they obtained be- 
tween the SFR and the stellar mass of the simulated galaxies 
(log(5 FR) = 0.66 log(M SMr )-5.69) gives SSFRs similar to those 
obtained with our model but with a steeper decrease toward high 
stellar masses. Finl ator et al.l d2006l) simulated B dropouts at z=4 
and also found a linear relation between the SFR and the stel- 
lar masses: log(SFR) = 1.141og(M SMr ) - 10.2. The resulting 
SSFRs are slightly lower than ours at z=4, and they also exhibit 
a flat distributi on of the SSFR w ith t he stellar mass . The r ela- 
tions found bv lDaddi etail d2007l) and lFinlator et all d2006l) are 
reported in Fig. [9] (panel b). 

6.3.2. Redshift evolution of the metallicity 

Chemical evolution models also allow predictions of the evolu- 
tion of the metallicity in star-forming galaxies. While it is not the 
case for our sample, recent surveys have allowed to determine 
abundances in relatively large numbers of high redshift galaxies. 
We thus think it is worth presenting our models predictions for 
the evolution of metals. 

In our models, we obtain a progressive increase in the metal- 
licity with time, at the same tie as the stellar mass increase with 
more massive galaxies always being more metal-rich (as a result 
of the star formation histories depending on the mass). At z=0, 
a metallicity-stellar mass relationship is present, corresponding 
to the we ll-known luminosity-met allicity relationship (e.g. the 
review bv lHenrv & Worthevill999l) . 

At higher redshift, we predict that a relation still exists (see 
Fig. [Tol l, but is progressively shifted to lower metallicities (by 
about 0.2 dex at redshift 1, 0.4 dex at redshift 2). We obtain in 
the models a steepening of the relationship between redshifts 1 
and 0: in the more massive galaxies, the metallicity increases 
more rapidly because the gas reservoir is no longer replenished 
by primordial infall, thus the newly formed metals are less di- 
luted even if the level of star formation is low. 

We note that the obtained change in metallicity is relatively 
modest as the theoretical uncertainties on the models are about 
a factor 2 (0.3 dex, but this is a systematic effect that should not 
affect our trends but only shift them, as a function of both magni- 
tude and redshift). Another uncertainty comes from the observa- 
tional side: it is indeed difficult to estimate abundan ces, and var- 
ious m ethods produce different results. For instance, Lian g et al] 
(2006) have shown that the equivalent width method produces 
systematically higher abundances by 0.2 dex. It is thus likely 
that part of the difference in the observed samples discussed be- 
low is due to the method used to measure O/H. These variations 
in metallicity should th en be addressed very carefully. 

iRupke et al.l d2007l) have recently found that the abundances 
in LIRGs increased by ~0.2 dex from z~0.6 to z~0.1, which is 



roughly consistent with our predictions. Nevertheless, we must 
remain cautious since the nature of the LIR Gs at low and in- 
terme diate redshift may well be different (e.g. Melb ourne et all 
120051) . 

There are several studies of the metallicities in high redshift 
galaxies, or of the luminosity-metallicity relationships evolution. 
However, due to various selection biases and various ways of 
estimating the metallicities, they are often difficult to compare 
with each other. For this re ason, no consensus has been found 
yet dLamareill e et all [20061 and reference therein). In this paper, 
they investigated the metallicities of 131 intermediate redshift 
star-forming galaxies. This allowed them to compare the local 
and intermediate redshift (0.2 <z< 1), split in 0.2 redshift bins) 
mass-metallicity relationship. They did not find any significant 
evolution of the slope, but did find that the high redshift relation- 
ship is shifted to lower metallicities (as predicted by our model) 
by 0.28 to 0.55 dex (depending on the analysis performed) at z 
~ 1 . This is only s l ightly more than the value we predicted. 

Sava glio et al.l d2005l) used the Gemini Deep Deep Survey 
and the Canada-France Redshift Survey to investigate the stellar 
mass-metallicity relationship between z ~ 0.1 and z ~ 0.7. They 
also found that the metallicity is lower at higher redshift for the 
same stellar mass by ~ 0.15 dex (according to their figure 13). 
They proposed a closed-box model for the chemical evolution 
of galaxies able to reproduce this result. This model is based on 
a simple exponential decline of the star formation rate after a 
formation redshift equal to 3. The dependence of the exponen- 
tial folding time on the total mass of the galaxy is fine-tuned to 
reproduce the shift they found in the mass-metallicity relation- 
ship (also taking the metallicities of distant Lyman break galax- 
ies into account). Such a model is similar in spirit to the one we 
propose, but ours is based on a more detailed modelization of 
the physics of galaxy evolution (infall, star formation law, finite 
lifetimes of stars) and was not fine-tuned to reproduce the evolu- 
tion of a simple property between low and high redshifts but of 
ma ny properties of nearb y gal axies. 

iHammeretal] d2005l) and iLiang et all d2006) also find that 
Z ~ 0.7 emission line galaxies were poorer in metals than 
present-day spirals, by 0.3 dex, a resul t also consistent with 
the a verage trend found by iMaier et al.l d2005l) . On the other 
hand, Kobulnicky & Kewlev (2004) report a smaller variation in 
0.14 dex from z = to z = 1. The factor two difference with 
other studies might be linked to the fact that they used equiv- 
alent widths, and standard underlying stellar absorption rather 
than using high quality calibrated spectra and measuring the 
Balmer absorption. Their data suggest a steepening of the slope 
of the metallicity-uminosity relationship. One should also note 
that a direct comparison in this case is difficult since the models 
show the metallicity-stellar mass relationship. While luminosity 
scales with the mass in nearby normal galaxies, the luminosity of 
high redshift galaxies might be str ongly affected b y their current 
star formation rate. Going further. Erb et al.l d2006l) suggests that 
star-forming galaxies at redshift ~ 2 have 0.3 dex fainter metal- 
licities. This is a somewhat smaller decrease than expected from 
our models; but taking into account the issues of selection and 
metallicity calibration, it is (at least) not inconsistent with them. 

Because of the uncertainties in measuring metallicities and 
in converting stellar masses to B band luminosities, a more de- 
tailed comparison is quite pointless. We hope that Fig. [TOl gives 
a sketch of the evolutionary trends expected in the case of sec- 
ular evolution, from whi ch eventual departu res should be due to 
more complex physics (Ru pke et all [2007). Determinations at 
high redshift of star formation rates, stellar masses, and metallic- 
ities may be compared to these trends to know how much galax- 



V. Buat et al.: Star formation history of galaxies from z=0 to z=0.7 



13 



8 



-9 



— a — 



E-io 

w 
o 



-11 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 11/ 




- 






- 'Sr 








- s ' / 




-// 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 



12 3 4 
redshift 



CO 

o 



-8 



-9 



10 



-b- 



-11 



I I I I Ml I I | I I I I | I I I I 

^ F06 , z=4 

KJ3"" 




I i i i i I i i i i I 



8 9 10 11 12 
log(M star ) (solar units) 



Fig. 9. Predictions of the models for higher z. All quantities are in solar units. Panel a: SSFR as a function of z for three rotational 
velocities and therefore masses (cf Table 1). Dotted line: 150 km s _1 , dashed line: 220 km s _1 , solid line: 360 km s _1 . Panel b: SSFR 
as a function of the stellar mass for different redshifts, each symbol corresponds to a specific redshift quoted on the plot, the values 
corresponding to a same redshift are connec ted with a red solid line. Dotted lines corresponds to the regression lines proposed for 
the SAMs of Finlator et all (12006) (F06) and lDaddi et al1 (l2007h CD07). 
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Fig. 10. Evolution of the metallicity as a function of z (panel a) and of the stellar mass (panel d), the symbols and lines are the same 
as in Fig. [9] 



ies are, on average, far from a secular evolution or not at each 
redshift. For instance, the values of st ellar masses and m etallici- 
ties in z — 5 Lyman break galaxies of lAndo et all (120071) are in- 
compatible with simple secular evolution objects. It is thus likely 
that these objects are not precursors of normal nearby spirals. 



7. Conclusions 

We have analyzed the star formation rate and the stellar mass of 
galaxies selected to be active in star formation at z=0 and z=0.7. 
The selection was performed in UV and in IR (rest frame). As 
long as relatively massive systems are studied (M star > 10 10 M G ), 
the IR selection is found to be more efficient than the UV one 



to select all the star-forming systems. Nevertheless the galaxies 
selected in UV and in IR exhibit similar variations in the SSFR. 

We compared mean relationships between the observed spe- 
cific star formation rate and the stellar mass at z=0 and z=0.7 
with physically motivated models aimed at reproducing the 
mean properties of local spiral galaxies and of the Milky Way. 
These models are based on a progressive infall of gas into the 
galactic disk starting at high z. The agreement is found to be 
quite good given the uncertainties in the models and the data. 
Both data and models exhibit a fairly flat distribution of SSFR 
for galaxies with masses between 1O 1O M and 1O U M , this flat- 
tening being more pronounced at z=0.7 than at z=0. These re- 
sults are consistent with those obtained at intermediate redshift 
from surveys selecting star-forming galaxies. 
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We have proposed an analytical formula for the mean evo- 
lution of the SFR with time and the parameters vary with the 
current stellar mass of the galaxy. This formulation will allow 
anyone to perform very simple calculations predicting e.g. star 
formation rates, stellar masses, or metallicities. We present pre- 
dictions for the values of the specific star formation rate, stellar 
mass, and metallicity at high redshift. These predictions can be 
taken as templates for a secular evolution of galaxies, including 
gas accretion. They are found consistent with the mean trends 
deduced from the simulations of semi-analytical models up to 
z=4. We tentatively compared these predictions to some existing 
data concerning star-forming galaxies up to z~2 without finding 
major departures. 

The comparison of model predictions with observations is 
difficult because of uncertainties (in observations, models, star 
formation rate, or metallicity calibrations) and selection biases. 
Thus we wish to stress that the trends presented in Fig [9] are 
only indications of the expected mean evolution in a simple sce- 
nario. Caution should be taken when comparing them to data. 
However, we believe they will still be useful for comparisons 
with the results of future or on going large surveys, especially 
if star formation rates, stellar masses, and metallicities are com- 
puted in consistent ways at different redshifts. 
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